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Abstract 

We present the program Top++ for the numerical evaluation of the total inclusive cross-section for producing 
top quark pairs at hadron colliders. The program calculates the cross-section in a) fixed order approach with 
exact NNLO for qq ^ tt + X, qq tt + X, qq' ti+ X, qq' ti + X and through approximate NNLO 
for gg ^ tt + X and b) by including soft-gluon resummation for the hadronic cross-section in Mellin space 
with full next-to-next-to-leading logarithmic accuracy. The program offers the user significant flexibility 
through the large number (31) of available options. Top++ is written in C-| — h. It has a very simple to 
use interface that is intuitive and directly reflects the physics. The running of the program requires no 
programing experience from the user. 



Program summary 

Name of the program: 

Program's homepage: 
License, Warranty: 
Compiler: 
Operating system: 
Program language: 
Memory required to execute: 
External libraries: 

Keywords: 

Typical running time: 
Accuracy: 



Top++ (ver. 1.3). 



http : //www . alexEinderinitov . com/software 



GNU Public License. No warranty given or implied. 

Developed and tested with GNU Compiler Collection's C++ compiler. 

Linux; Mac OS X; can be adapted for other unix systems. 

C++. 

Typically less than 200 MB. 

GNU Scientific Library (GSL); the Les Houches Accord PDF Interface 
(LHAPDF). 

Top-quark, Resummation, QCD, Precision Physics, Hadron Colliers. 

Depending on the options. The program is optimized for speed. 

Sub per-mill accuracy achievable in realistic time (program does not 
employ Monte Carlo methods). 



♦Preprint numbers: CERN-PH-TH/2011-303, TTK-11-58 
Preprint submitted to Elsevier 



July 4, 2012 



What is new in ver.1.3: changes from ver.1.2. 

• Added the complete NNLO results [ij for the reactions: 
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Vqq. 



Adding the reaction qq ^ tt + (Z^jj^g (as defined in Ref J_^) completes the calculation of the reaction 
qq^ tt + X aX NNLO. Please note that, compared to version 1.2, version 1.3 adds few seconds to each 
evaluation due to the initialization of the newly added partonic fluxes qq, qq' ,qq' . 

• Added an option that allows switching between ".LHgrid" and ".LHpdf ' files in the LHAPDF interface. 
The option is called PdfFileType with default value LHgrid. The alternative value is LHpdf. 

• Added the symmetric prescription for computing PDF uncertainties (defined as in Ref [2]). It is 
automatically invoked for the ABMll [2], ABKM09 d and earlier Alckhin [i^ family of PDF sets. 



What is new in ver.1.2: changes from ver.1.1. 

With respect to version 1 . 1 (minor bug fixes with respect to version 1 . 0), Top++l . 2 includes the exact 
NNLO result for qq ^ tt + X. This is a major upgrade to the program's functionality that necessitates the 
following modifications: 

• Remove the option Cbarqq which is now obsolete due to the inclusion of the exact NNLO qq ^ tt + X 
result. 

• The hmctionalities of the following options are modified with respect to ver.1.1 and lower: DrderFO, 
OrderRES, RESonORoff, H2qq, H2ggl, H2gg8, NNLO, NNLOonORoff , Cbargg, ECMLHC. 

• The description of the option MtopLimit is improved. 

• Included is a new set of examples, that reflects the new functionality. 

• Slight modifications to the Manual that reflect these changes. 

• Added functionality to the file top++ . cpp that restricts the option MtopStep to be positive (see its 
description below), and prevents the program from possibly entering into an infinite loop. 

The default settings of the program reproduce the Tevatron NNLO numbers of Ref. The best 
recommended setting for the LHC can be found in the file top++best-LHC8 . cf g. They are a combination of 
the findings of Ref [a] with the exact NNLO result for the qq ^ tt + X reaction from Ref Q . For the sake 
of documenting the results, in the following we explicitly present our best predictions for the Tevatron, 
LHC @ 7 TeV and LHC @ 8 TeV: 

^NNLO+NNLL(Tevatron; = 173.3 GeV) = 7.067 '^^^'^^ [scales] tlf.l'S^, [pdf] , (5) 

^NNLO+NNLL(LHC,Tcv; m, = 173.3 GeV) = 156.2 t'.^to'S' t^^^M -t^JfrS [pdf] , (6) 

^NNLO+NNLL(LHCsTeV; = 173.3 GcV) = 222.5 tllZiri -5:6r2:5%) [P^f] ■ 



^We take this opportunity to update the legacy NLO+NLL numbers given in Tables 2,3 of Ref. 0]. The fully updated 
results, which can be obtained with any version of Top++, read: 172.0 i^ig* (7;9%] (Table 2, line 2); 165.8 ^i4;i|8;5%) (Table 2, 

line 3); 7.039 422 (6 o%) ('^^'^^^ ^' 6.911 ^0 497(7 2%) (Table 1, line 3). We thank Christian Schwinn for bringing the 

need for this update to our attention. 
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The numbers for the Tevatron Eq. ^ reproduce the best prediction of Ref. Q , while the ones for the LHC 
Eqs. (j6l7p are new, and represent the current state-of-the-art. 

The dramatic improvement from the inclusion of the NNLO corrections to qq —¥ tt + X at the Tevatron 
has already been discussed in Ref. Q . The inclusion of the NNLO correction in qq tt + X improves the 
theoretical prediction at the LHC with respect to the NLO-^NNLL one for both LHC @ 7 TeV and 
LHC @ 8 TeV the scale uncertainty bands perfectly overlap, and the total scale uncertainty in Eqs. (|6l7p is 
reduced, in absolute size, by about 2%. 



What is new in ver.1.1: changes from ver.1.0. 

With respect to version 1 . 0, Top++l . 1 fixes two bugs in the file Resummation. cpp. 

• Corrects the normalization of the hard matching constants H2qq, H2ggl, H2gg8 (the factor of I/tt^ 
was missing in ver.1.0). 

• Fixes a bug in the calculation of the matching of the all-order soft exponent. The bug was affecting 
the results in the unusual cases when lower logarithmic accuracy was requested than the fixed order 
one, i.e. the following three cases: NLO-^LL, NNLO-I-LL and NNLO-hNLL. 

These bug fixes do not affect any of the numbers published in this manual or in Ref. 6] . 



1. Introduction: what is this program for? 

This program calculates the top-pair total inclusive cross-section in hadron collisions. It has built-in 
the most advanced theoretical results available in the literature. The program is able to calculate the total 
top-pair cross-section in both pure fixed order perturbation theory and by including soft-gluon resummation 
through next-to-next-to- leading logarithmic order (NNLL). Top++ is the first publicly available program that 
can perform soft gluon resummation in top-pair production. 

The program is written in C-I--I- in a modern, modular and object-oriented way. It should be very easy 
to install on most Linux systems; please consult Appendix A| for details. Moreover, once installed, the 



program is trivial to run. The program has been written with a user in mind that has no programming 
experience whatsoever. For that reason, the program has a very simple user interface that is the only point 
of contact between the program and the user. The user interface is described in Section [21 Users with 
average programming experience will find the program very easy to customize to their own needs. That 
may not be needed, however, since all options that have arisen in our own work on the subject have already 
been pre-programmed and are easy to access directly through the user interface. 

In this manual we do not describe the physics in detail. A short description, needed to make the reading 
of the program's options self-consistent, can be found in Section HI The relevant description of the physics 
and the options implemented in this program can be found in Refs. 0, The results in these papers 
represent the most advanced results in top physics to date. Given the fast pace of developments in the 
subject, however, we expect that the complete set of NNLO results will be available soon. 

And an important word of caution regarding the interpretation of numerical results derived with the help 
of the program. Clearly a program with a number of options (like ours) is capable of correctly producing 
almost any number. To be able to interpret such a number, however, one needs to know exactly how it was 
derived. This can become a problem with the state of top physics where the most advanced results used 
to constitute a series of approximations (see Ref. @) for both Tevatron and LHC. To reduce the chance of 
mis-communicating results derived with this program, we have implemented default settings that correspond 
to the best prediction of Ref. Q for the Tevatron and a combination of Refs. 0, Q for the LHC, and we urge 
the users to use them for publications. Parameters that are not considered part of these settings are the 
parton distribution functions (PDF) and the value of the top quark mass. These should always be quoted 
to ensure complete documentation of the results. Other settings, mainly related to the workings of the 
program, can affect the final numerical value (through numerical precision) but that should be at a precision 
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level below what is physically relevant. As a rule of thumb we urge the user to produce numbers that are 
correct at the level of one per-mill (i.e. 0(10"'^)). Precision better than that is a pure matter of taste on 
the side of the user (and, at a certain level, of the ability of the otherwise very capable integration routines). 

2. How to use the program 

In the following we assume that the program has been installed and is running correctly; see [Appendix A 
for more details. Next, we assume that the program is installed in a directory called top + + and we are 
already there, i.e. the command line reads: □ 

~ top++$ 

Next one needs to open the file top++.cfg in a text editor. Using, for example, pico one types: 
top++$ pico top++.cfg 

Once the user has set all options at their desired values, the file needs to be saved. A number of examples 
are supplied with the program; see [Appendix B[ Then the user needs to execute the program by typing: 

~ top++$ ./top++ 

The program starts running, displays its step-by-step progress, timing pre-defined milestones. Once the 
run is completed, the final result is displayed and the program exits. The output can be found on the screen 
and in the file top-h+.res located in the program's directory. If a new run is desired, one has to simply 
repeat all the steps described above. 

In the following we describe all 30 options that are available to the user through the file top++.cfg. The 
options are grouped into five subgroups. Please note that since all options have predefined default values, 
the user needs to only specify the values for options that differ from their default values. In particular, if all 
options are at their default values the file top++.cfg can be empty. It might be convenient to keep certain 
options typed in the file, but when not in use the user can comment them out by putting the symbol / at 
the beginning of a line. 

1. General Setup (type of collider, PDF set, pure fixed order calculation versus one with resummation). 

(a) Collider: Takes two values: TEV (default) or LHC as labels for pp or pp colliders. 

(b) WithResummation: If set to YES (default), the program will compute the observable at N"LO + 
N^LL depending on the values for n and k specified in the subgroup Resummation below. In 
this case all settings in the subgroup Fixed Order become irrelevant. 

(c) PDFuncertainty : If set to YES the program will compute and display the PDF uncertainty. 
Three methods for computing PDF uncertainties are available: First, the symmetric prescription 
(as defined in Q) is automatically invoked for the ABMll [ij, ABKM09 Q and earlier Alekhin 
[3] family of pdf sets. Second, the prescription for the NNPDF family of sets 0] is selected auto- 
matically for any NNPDF set and, third, for all other PDF sets the program uses by default the 
asymmetric prescription of Ref. Q (see also Ref. [3). In particular, the latter prescription is the 
most appropriate one for the MSTW2008 family of sets jlO| . When the PDF uncertainty is being 
calculated, central scales choice fip — fJ-R — "nit is set automatically, independently of the values 
for these scales set by the user (as described below). If PDFuncertainty is set to NO (default) 
then the program computes and displays the scale variation for the choice of fiF,fiji specified by 
the user (explained below) for a single pdf member (specified by the user; see option [TJ below). 
Please note that the definition of a "central" set depends on the PDF set. 



^Please note that the precise text on the left of the cursor depends on the particular terminal. 
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(d) RestrictedScaleVariation: Takes either NO or a number greater than or equal to 1.0. This 
option allows the user to set any restriction (or no restriction at all if RestrictedScaleVariation 
NO) on the allowed ratio for the renormalization and factorization scales. For more information see 
the description of the group of options Top quark mass and renormalization/factorization 
scales. The default value is 2 . which corresponds to the restricted scale variation of Refs. [5|,Q. 

(e) PDFset: The PDF set. The program uses the LHAPDF library and follows its nomenclature. 
The default set is MSTW2008nnlo68cl. Both ".LHgrid" and ".LHpdf files can be used, see the 
option Pdf FileType below. 

(f ) PDFmember : The specific member of the PDF set PDFset the user wishes to use for the calcu- 
lation of the cross-section. Please note that the definition of a "central" set depends on the PDF 
set. Default value is 0. 

2. Top quark mass and renormalization/factorization scales Please note that the ranges of muR 
and muF (defined in the following) that the user may provide need not be equal in length. The 
range of their ratio is unrestricted if RestrictedScaleVariation NO is chosen. When the option 
RestrictedScaleVariation takes a numerical value (with the constraint RestrictedScaleVariation > 
1.0), then the ratio of the two scales is restricted between: 

1 

^ I^f/i^r < RestrictedScaleVariation. 

RestrictedScaleVariation 



(a) Mtop: The value of the on-shell top mass (in GeV). If the user would like to loop over a range 
of values for the top mass (see the description of the following two options) Mtop represents the 
lower end of that range. Default value is 173.3. 

(b) MtopLimit : The upper limit of the range of values for the top mass that the user wants to loop 
over. If MtopLimit < Mtop then the program automatically sets MtopLimit = Mtop and issues 
a warning message. The default is the unattainable value MtopLimit = —1. The value —1 is 
converted internally to MtopLimit — Mtop too, but is special, because only for it the (annoying) 
warning message is suppressed. 

(c) MtopStep: A positive number (need not be an integer) that specifies the step with which the 
value of the top mass is incremented in a loop. Its default value is 1. 

(d) muR: Set of values for the renormalization scale (in units of the top mass Mtop). By default it 
takes the value 1.0. It can take an arbitrary number of values. In particular, to determine the 
scale variation of a cross-section, one may use the set: 

muR 0.5 1.0 2.0 . 

A fine scan can be achieved with the help of the following set: 

muR 0.50 0.55 0.60 0.65 0.70 0.75 0.80 0.85 0.90 0.95 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.0 . 

The larger set takes much longer to compute and typically returns the same result as the restricted 
set consisting of three elements. We have encountered, however, exceptional situations that in 
our experience lead to differences of up to 0.5%. 

(e) muF: Set of values for the factorization scale (in units of the top mass Mtop). It is set and used 
independently of the renormalization scale described above. Its default value is 1.0 and can be 
set analogously to muR. 

3. Resummation (all options in this group are irrelevant when WithResummation NO is chosen). 
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(a) OrderFO: Takes values LD, NLO or NNLD (default). It specifies the order n of the fixed order 
cross-section that should be added to the resummed result, i.e. N^LO + N^LL. For example LO 
implies n = 0, NLO implies n = 1, etc. The option NNLO includes the exactly known NNLO result 



for qq ^ tt + X and the approximately known NNLO cross-section for gg ^ tt + X; see also 
the group of options Fixed Order below. 

(b) OrderRES: Takes the values LL, NLL and NNLL (default). It represents the logarithmic order k 
of the resummation, i.e. N'^LO -I- N^LL. For example LL implies fc = 0, NLL implies fc = 1, etc. 
Full NNLL resummation is implemented matched to a hard function with H2qq — 84.81 and 
the still unknown coefficients H2ggl , H2gg8 set to zero by default (see also the description of the 
options below). 

(c) A: The value of the parameter A introduced in Ref. [l2]. Typically used with the values A 2 
(default) and A 0. 

(d) RESonORoff : Takes the values ON or OFF (default); specifies how the unknown two- loop hard 
function coefficient H2ggl , H2gg8 are treated. It does not affect H2qq. Wc have implemented this 
option consistently both for the fixed order (relevant if OrderFO is set to NNLO) and the resummed 
components. Moreover, the unknown constant Cgg in the fixed order NNLO approximation is 
automatically expressed through H2ggl , H2gg8 (described in option [Sj- 

(e) TwoLoopCoulombs : Takes the values YES (default) and NO. It includes/excludes the two loop 
Coulombic terms in the function 

^(Coui)^ see Eq. ® below. 

(f) H2qq, H2ggl, H2gg8: The three constants in the two-loop hard function. They are defined 
in a normalization a^/Tr. H2ggl, H2gg8 are irrelevant when RESonORoff is set to OFF. When 
RESonORoff is set to ON, one needs to vary them to determine the full theoretical uncertainty. 
The default values for H2ggl and H2gg8 are 0.0 and for H2qq is 84.81 

4. Fixed Order (all options in this group are irrelevant when WithResummation YES is chosen). 

(a) LO: takes values YES (default) or NO. 

(b) NLO: takes values YES (default) or NO. For example, to perform a calculation in NLO QCD one 
needs to set LO YES and NLO YES. This way the user has direct access to the individual orders in 
the perturbative expansion. Same applies to the option NNLO below. 

(c) NNLO: takes values YES (default) or NO. This option implements the exact NNLO result for 



qq ^ tt + X and the threshold approximation to the NNLO cross-section for gg ^ tt + X 



(d) NNLOonORoff : Takes the values ON or OFF (default). Specifies how the ln(/i) terms in the 
approximate NNLO result in gg ^ tt + X are treated; see Ref. Q. Irrelevant for qq ^ tt + X. 

(e) Cbargg : The value for the yet- unknown two-loop constant Cgg in the threshold expansion of the 
NNLO gg ^ tt + X cross-section [ll| (defined in a normalization as/(47r)). The default value 
for this constant, to be used at the Tevatron, is Cbargg = 0.0. As explained in Ref. @ a more 
appropriate value for the LHC is Cbargg = —1136.81 (corresponding to H2ggl = H2gg8 = 0). 

Setup parameters (parameters related to the working of the program and other, less frequently 
modifiable parameters). 

(a) ECMLHC: The cm. energy of the pp cohider (in GeV). The default value is 8000. 
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(b) 



ECMTEV: The cm. energy of the pp colhder (in GeV). The defauh value is 1960. 



(c) 



Precision: Defines the required relative precision of the integration routines: 



Relative precision =10 



Precision 



The default value precision 2 tends to produce fast and accurate (at the per-mill level) results. 

(d) NPdfGrid: Defines the size of the grid on which the PDF fiuxes are being discretized. The 
default value NPdfGrid 100 tends to produce fast and sufficiently accurate results. 

(e) ETA : Parameter that controls the subtraction flux implemented in the resummed calculation. 
The default value ETA le-5 is optimal. H The user normally will not need to be concerned with 
this parameter. 

(f) CMP: The value where the contour for the Mellin inverse transform crosses the real line. The 
default value CMP 2 . 7 is optimal. The user normally need not be concerned with this parameter. 

(g) PdfFileType: Option for switching between ".LHgrid" and ".LHpdf flies in the LHAPDF 
interface. The default value is LHgrid. The alternative value is LHpdf. 

3. Once the program is running: some fine tuning 

3.1. Numerical precision and speed 

As every program for numerical calculations, Top++ has its limits, too. In the following we discuss this, 
as well as ways to improve the accuracy and length of the runs. 

There are two places where speed (and therefore accuracy) can be controlled. The first one is through 
the option Precision of the integration routines. In our own experience the value Precision 2 is more 
than adequate to calculate the tt cross-section to per-mill accuracy. Increasing the value of Precision slows 
the calculation down. The program is set in such a way, that if the integration routines cannot reach the 
accuracy in a point, a warning message is displayed specifying the relative error returned by the integration 
routines. The presence of these messages, as such, is harmless. They can be ignored if the displayed relative 
precision is high enough. Rarely, the output might contain the symbol NaN, which stands for Not-a-Number. 
NaNs are returned due to invalid numerical operations such as division by zero, or operations, which produce 
numbers out of the range of double precision. This is a sign of numerical instability. We have taken the 
pragmatic approach of not handling such exceptions in any specific way (which is very difficult in practice). 
In case of a NaN result, the user should rerun the calculation with increased requested accuracy. 

We recommend that at least once the user does the calculation with a larger value for Precision and 
verifies that the change in the result is beyond the required accuracy. 

A second (and independent) source of numerical uncertainty is the size NPdfGrid of the grid over which 
the partonic fluxes are approximated. We have implemented a second order flnite difference scheme. The 
relative precision scales as ^ 1/NPdf Grid-^. In practice we have noticed that a value NPdfGrid 100 is more 
than adequate in terms of accuracy and produces very fast calculations. We recommend that at least once 
the user does the calculation with a larger value, say NPdfGrid 500, and verifles that the change in the 
result is beyond the required accuracy. 

Finally, the user should keep in mind that the overall numerical accuracy is a combination of the settings 
for Precision and NPdfGrid, i.e. increasing only one of them to an extreme may not have a net positive 
effect on the overall uncertainty but might lead to a significant slowdown. 



^We use the usual notation: le-p = 10 
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3.2. Outputting the results 

The program outputs on the screen the information about the timing of each step and aU results. The 
specifics depend on the requested options. At the end of the calculation, a summary of the final result is 
displayed, including the central value (i.e. the value corresponding to central scale choice, if requested) and 
the scale or PDF variation's absolute and relative values. The program also outputs the final result in a file 
top++.res which is ready for plotting with the program Gnuplot. This is particularly useful for the case 
when the user requests a loop over a range of values of the top quark mass. In such case the result for each 
value is conveniently recorded in the file. Please note that the file top++.res is overwritten after each run. 

3.3. Modifying the default values of the parameters 

The default values for all 30 parameters available to the user are set inside the function mainO located 
in the file top++.cpp. All default settings can be modified by the user, although this is not recommended 
and should not be necessary. 

3.4- Additional 

The strong coupling constant is calculated at a scale through the LHAPDF interface The 
calculation of the cross-section is performed in a scheme with Np = 5 active fiavors. For consistency only 
PDF sets with Np = 5 active flavors should be used for scales above nitop. 



4. Contact with physics 

The program Top++ comes with settings that reproduce the "best" numbers from Ref. [5| for the Tevatron. 
The best numbers for the LHC will now slightly differ from the ones in Ref. Q due to the improved precision 
after the inclusion of the NNLO result for qq tt+X. These are the best current theoretical predictions with 
realistic theoretical uncertainties. All physics-related inputs and approximations are defined and explained 
in Ref. [6], which should be consulted if more detailed questions arise. 

The program calculates the tt total inclusive cross-section, schematically given as: 

^tot = 4"o. + - ^"lag] ■ (8) 

The power n is set by the option OrderFD if resummation is desired (i.e. when WithResummation YES). 
The convention is that n — 2 corresponds to LO, n = 3 to NLO and n — 4 corresponds to NNLO. If 
WithResummation NO, i.e. one wishes to calculate in a fixed order approach, then the terms in the square 
bracket in Eq. ([5]) are absent and n is set through the options LO, NLO and NNLO. 
In Mellin TV-space, the resummed partonic cross-section reads 

res part \ ^ (Coul) (Hard) a 

1=1,8 

The function cr^^°"'^ contains the Coulombic effects and has a known perturbative expansion through 
NNLO. The function CT(Hard) jg 7V_independent hard function. The NNLO constant terms H2ggl , H2gg8 
are currently unknown. These constants, in an expansion in as /t^, correspond to the options H2ggl , H2gg8. 
The option H2qq is available to the user, too. However, since it is set by default to the exact value derived in 
Ref. [5[, normally it should not be modified. The ln(/i) terms in the two- loop 5(7— contribution to the function 

^(Hard) g^j.g 

included if RESonORof f ON and excluded when RESonORof f OFF. The latter option gives a more 
conservative (i.e. realistic) determination of the theoretical uncertainties through scale variation alone. 
The function A^v.i contains the towers of LL, NLL and NNLL soft ln(A'') logs. 

By setting the option OrderRES LL the program included only the LL corrections from A^rj, as well as 
the LO terms in the functions cr(Coui) ^^^^ ^(Hard)^ OrderRES NLL then Aatj includes LL and NLL logs, 
while the functions (t(*^°"') and a^^^^'^^ are given by their NLO expansions. Similarly, for OrderRES NNLL 
Ajvj includes LL, NLL and NNLL logs, and the functions (t(^°"') and a^^''"^^ are given by their expansions 
through NNLO. The expressions for these functions can be found in Refs. 
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For pure fixed order calculations one also has a variety of options for the treatment of the currently 
approximately known NNLO corrections in gg ti + X. Setting NNLOonDRof f ON one includes the two- 
loop ln(/i) terms unaccompanied by a power of /? or ln/3. By setting NNLOonORoff OFF one excludes these 
two-loop logarithmic terms which corresponds to the standard setting adopted in Ref. Q. Note that the 
qq ^ tt + X reaction is unaffected by this option. 

(2) 

Through the options Cbargg one has a direct access to the currently unknown two- loop constant Cgg . 
This is helpful for determining the sensitivity of the approximate NNLO cross-section to this unknown 
parameter. 



5. Summary 

We present the C-1--I- program Top++ for the calculation of the total top-pair cross-section at hadron 
colliders. This is the first publicly available program capable of performing soft-gluon resummation for this 
collider observable. The program uses all currently available theoretical results. The user has access to 30 
options which results in a great deal of fiexibility and control over the calculation on the side of the user. In 
this manual we have given only a very short introduction to the physics behind our program. For further 
details the user should consult Ref. [6|. 

The program is organized in a modular, object oriented way. It can seamlessly accommodate future 
results, notably the remaining NNLO results, by combining them with the currently available NNLL resum- 
mation and known NNLO results for qq ^ tt + X 

The program is optimized for speed given the user's requirements for accuracy. Our experience shows 
that, depending on the chosen options, the run times can vary significantly. In practice the speed can be an 
issue only for resummed calculations, due to the integration of rapidly oscillating functions in the complex 
plane. For fixed order calculations the run times are very short. 

As we already emphasized Top++ is the only publicly available program capable of producing resummed 
calculations. Another program for the calculation of the total cross-section in fixed order perturbation 



theory is available in the literature While we have not performed detailed side- by-side comparisons, 
our experience shows that our program can be significantly faster, especially when higher accuracy is required 
(since we do not use Monte Carlo integration methods) . We would like to also note that a direct comparison 
between the results of Top++ and HATHOR ^14*1 for the approximate NNLO results may not be possible. 
The reason is that different implementations of this approximation are contained in both programs. In 
Top++ we have implemented two versions of this approximation, depending on whether the two-loop scale 
dependent terms are included or excluded (see Ref. [gI for details). This is controlled through the option 
NNLOonORoff OFF/ON. An approximation analogous to our NNLOonORoff ON is implemented in HATHOR by 
using the option LOG_ONLY. We note that while the two are compatible they are not identical, since different 
subleading terms 0{f3) are included in both implementations. 

We would also like to stress (see Ref. @) that approximations like NNLOonORoff ON could significantly 
underestimate the theoretical uncertainty if its estimate is based on scale variation alone. 
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Appendix A. Installation 

The program has been written in standard C++ and has been tested to correctly compile under the GNU 
compiler g++ from version 4 upwards. It requires two external libraries: 
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• GNU Scientific library, which can be downloaded free of charge from 
http : / / www . gnu . org/s/gsl/ 

and is used for special functions and integration 

• the Les Houches Accord PDF Interface, which can be downloaded free of charge from 
|http : //pro j ects . hepf orge . org/lhapdf /| 

and is used for the parton distribution functions 

To setup the program for installation, it is necessary to set three variables in the Makefile contained in 
the Top++ installation directory. 

• CXX - C++ compiler, specify the full path if necessary 

• GSLDIR - prefix directory for the gsl library, it is assumed that the library is in GSLDIR/lib and the 
include files in GSLDIR/include 

• LHADIR - prefix directory for the Les Houches PDF library, it is assumed that the library is in 
LHADIR/lib and the include files are in LHADIR/include/LHAPDF 

In the next step, it is sufficient to compile the code using 
top++$ make 

The program should compile without any error massages or warnings, and is ready to use under the name 
top++. 

Alternatively to modifying the Makefile it is possible to compile directly with 

~ top++$ make CXX="user CXX value" GSLDlR="user GSLDIR value" LHADlR="user LHADIR value" 

The quotation marks above are not necessary, unless the paths contain spaces or other special values (as 
usual under unix). 

Appendix B. Examples (updated for ver.1.3) 

A number of examples can be found in the directory examples. The user can copy/paste their content 
into the file top++.cfg located in the program's main directory. 

Example 1: If one executes the program with its default settings (provided in the file top++def ault . cf g) 
one should get the following result (this result can also be obtained with an empty file top++.cfg or, alter- 
natively, one with all its options commented out): 

sigma_tot = 7.06698 [pb]. 

Example 2: With the file top++best-TEV. cf g one can compute the scale variation of the cross-section 
corresponding to the best prediction of Ref. [5] for the Tevatron: 

sigma_tot = 7.06698 + 0.143258(2.02715%) - 0.231754(3.27939%) [pb]. 

Example 3: Wc recommend the setting in the file top++best-LHC8 . cf g for the computation of the 
cross-section for LHC @ 8 TeV (or any other LHC energy after adjusting the value of ECMLHC). It is based 
on the results of Refs. 0, Q and reads: 

sigma_tot 222.664 + 21.2071(9.52428%) - 11.3969(5.11845%) [pb]. 

This result perfectly overlaps with the LHC @ 8 TeV number quoted in Eq. (26) of Ref. @ and, as anticipated, 
has about 2% smaller total scale uncertainty. 
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Example 4: To derive the PDF variation for the best prediction of Refs. Q at the Tevatron one can 
use the file top++pdf var-TEV. cfg: 



sigma_tot = 7.06698 + 0.186404(2.63768%) - 0.122316(1.73081%) [pb]. 



Example 5: To derive the PDF variation for the best prediction of Refs. at the LHC @ 7 TeV one 
can use the file top++pdf var-LHC7 . cfg: 



sigma_tot = 156.272 + 4.12863(2.64196%) - 4.19779(2.68622%) [pb]. 



Example 6: Fixed order calculations can be done by turning off the option for resummed calculation. 
For example, the scale variation of the NLO result at the Tevatron, computed for consistency with an NLO 
PDF set (MSTW in this case), is obtained from the file top++NLO-tev. cf g: 



Example 7: The fixed order, partially NNLO result for the LHC @ 7 TeV, with excluded two-loop 
scale dependent terms in gg ^ tt + X Q and non- vanishing two-loop constant C^g , is obtained from the 
file top++Fixed-Drder-LHC7.cfg: 



This result overlaps very well with the corresponding number quoted in Table 2 of Ref. 6] and has, as can 
be anticipated, about 1% smaller total uncertainty due to the inclusion of the exact NNLO for qq — > ti+ X. 

(2) 

Example 8: The fixed order NNLO result at the Tevatron, with vanishing two-loop constant Cgg , is 
obtained from the file top++Fixed-Order-TEV. cfg: 

sigma_tot = 7.00564+ 0.202329(2.88808%) - 0.30846(4.40302%) [pb]. 

This result corresponds to Eq. (11) in Ref. 

Appendix C. Program's structure: a brief overview 

The function mainO is located in the file top++.cpp. The program consists of four classes that are 
initialized in the function mainO . The classes, listed in the order they are initialized, and their functionality 
are briefly described in the following. 

1. Class PartonicFlux . This class precomputes the partonic fluxes on a grid of NPdf Grid points. One 
object of the class is created for each value of the factorization scale hf- See also Section IXTl 

2. Class FixedOrder. This class represents the partonic fixed order cross-section. One object of this 
class is created for each combination of {^Ft^-r)- 

3. Class SubtrFlux. This class implements a fake partonic fiux that mirrors the actual partonic fiux 
within a distance ©(ETA) from the partonic threshold. Our implementation follows Ref. [l3| where 
this flux was introduced and the need for it explained. Objects of class SubtrFlux are created for each 
object of class PartonicFlux. 

4. Class Resummation. One object of this class is created for each pair (fiF,fJ-R) if a resummed calcu- 
lation is requested by the user. Each object is constructed analytically in Mellin iV-space and then 
inverted numerically back to a;-space with the help of the Minimal Prescription of Ref. [l^ ■ 

A small number of functions can be found outside the above classes. The file Igamma contains the logarithm 
of the Euler gamma function Inr(z) for complex argument, and psin contains the polygamma function 
^(z,fc), fc = 0, 1 for complex z. The remaining functions are located in the file Utilities. These are 
the prescriptions for calculation of PDF uncertainties (we have implemented two prescriptions; see the 
description of option PDFuncertainty in Section [5]) and the actual computation of the final result as a 
convolution of perturbative functions and partonic fiuxes. 



sigma_tot = 6.68123 + 0.35577(5.32492%) - 0.752327(11.2603%) [pb]. 



sigma_tot = 154.779 + 11.0044(7.10975%) - 8.59923(5.55581% 



) [pb]. 
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